SPECIFICATION 



TO ALL WHOM IT MAY CONCERN: 

BE IT KNOWN that we, Hiroo OHMORI , a subject of Japan 
and residing at Yokohama-shi, Kanagawa, Japan, Takahiro ASAI a subject of 
Japan and residing at Yokosuka-shi, Kanagawa, Japan and Tadashi 
MATSUMOTO, a subject of Japan and residing at Ebina-shi, Kanagawa, 
Japan have invented certain new and useful improvements in 

"ADAPTIVE EQUALIZATION METHOD AND ADAPTIVE 
EQUALIZER" 

and we do hereby declare that the following is a full, clear and exact 
description of the same; reference being had to tiie accompanying drawings 
and the numerals of reference marked thereon, which form a part of this 
specification. 
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ADAPTIVE EQUALIZATION METHOD AND ADAPTIVE EQUALIZER 

BACKGROUND OF THE INVENTION 

The present invention relates to an adaptive equalization method and 
5 an adaptive equalizer for use, for example, in mobile communications to 
reduce the influence of fading. 

In a mobile communications environment signals transmitted from 
mobile stations arrive at a base station from various directions under the 
influence of surrounding buildings, trees and so forth. The influence of 

10 intersymbol interference, which results from the arrival of paths of large delay 
time differences in such a multipath channel, poses a serious problem for 
mobile communications. A known technique for lessening the influence of a 
signal waveform distortion by the intersymbol interference is an equalization 
technique that estimates the characteristic of the channel and performs 

1 5 equalization of delayed paths based on the estimated characteristic. One of 
equalization algorithms used in the equalization technique is a 
turbo-coding-technique-applied iterative equalization algorithm that achieves 
a very excellent error rate characteristic in low SNR. The turbo coding 
method implements, with practical coding-decoding complexity, an error rate 

20 characteristic that is almost a ceiling value obtainable with a channel coding 
method; only the presence of such an error rate is merely suggested in 
information theory. With the turbo coding method, the transmitting side 
transmits coded information sequence and its rearranged or interleaved 
version, and the receiving side repeatedly decodes them alternately with each 

25 other. Since the process of generating a coded sequence of convolution 
codes in the turbo coding is equivalent to the process of generation of the 
intersymbol interference in multipath channel, the decoding process of the 
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turbo coding technique can be applied to tihe equalization of delayed paths. 

Fig. 1 illustrates in block form an example of a transmitter-receiver for 
the iterative equalization disclosed in Daryl Reynolds and Xiaodong Wang, 
"Low Complexity Turbo-Equalization for Diversity Channels," 
5 http://ee.tamu.edu/~ reynolds/ . A transmitter 10 encodes an information bit 
sequence d(i) by an encoder 11, then rearranges the coded output sequence by 
an interleaver 12, and modulates the carrier by a modulator 13, thereafter 
transmitting the modulated output bit sequence b(i) over a channel 14. In a 
receiver 20 delayed paths are equalized by a SISO (Soft-Input Soft-Output) 
!i 10 equalizer 21 for the received signal r(t), and the log-likelihood ratio (LLR) Ai 
[3 between +1 and -1 for each code bit b(i) of the received signal r(t) is derived 

' J by the following equation (1 ). 

m 

fi A^[b(i)] = log™^±lll(t)l 

''^^^ ^Pr[b(i) = -l|r(t)] 

!U ^ p[r(t)|b(i) = +l] ^ Pr[b(i) = -H] 

^p[r(t)|b(i) = -l] ^Pr[b(i) = -1] 

i =:^i[b(i)] + X2[b(i)] (1) 

In Eq. (1) the first term, ?^i[b(i)], is extrinsic or external information 
15 derived by subtracting a priori information A-2[b(i)], by a subtracter 26, from 
the log-likelihood ratio Ai[b(i)] provided from the SISO equalizer 2 1 . The 
second term, A-2[b(i)], is the a priori information that is a rearranged or 
interleaved version, by an interleaver 25, of A,2[b(n)] obtained by subtracting 
A.i[b(n)], by a subtracter 24, from a log-likelihood ratio A2[b(n)] provided 
20 from a SISO decoder 23 in the previous iteration. The initial value of 

A-2[b(n)] is zero. The a priori information A-2[b(i)] is used as an a priori LLR 
value of a coded bit string {b(i)} . The extrinsic information is restored by a 
de-interleaver 22 to ?ii[b(i)], which is input as a posteriori information to the 



-3- 

SISO decoder 23 . The SISO decoder 23 uses the a posteriori information to 
derive a a posteriori LLR value A2 of the coded bit string {b(i)} by the 
following equation (2). 

A2[b(n)]=X2[b(n)]+:Vi[b(n)] (2) 
5 In the second and subsequent iterations, X,i[b(n)] is subtracted by the 

subtractor 24 from the a posteriori LLR value A2 derived by the SISO decoder 
23 to obtain extrinsic information ?.2[b(n)], which is provided as a priori 
information 5\.2[b(i)]via an interleaver 25 to the SISO equalizer 2 1 . Such 
iterations of equalization and decoding update the LLR value, providing a 

10 more probable LLR value. That is, the first log-likelihood ratio Ai is 

calculated by the first adaptive equalization of the received digital signal, and 
the first log-likelihood ratio Ai is decoded to obtain the second log-likelihood 
ratio A2, which is used to perform the second adaptive equalization of the 
received digital signal to obtain the first log-likelihood ratio Ai, which is then 

1 5 decoded, and the decoded value is output. 

In the iterative equalization the SISO equalizer and the SISO decoder 
use a MAP (Maximum A-posteriori Probability) algorithm, Log-MAP 
algorithm based on the MAP algorithm, or Max-Log-MAP algorithm that has 
reduced the computational complexity of equalization and decoding through 

20 use of approximate calculation. Still another algorithm is SOVA 

(Soft-Output Viterbi Algorithm) that provides bit reliability information while 
performing optimum decoding on a sequence-wise basis. In the case of 
using any of these algorithms for equalization by the equalizer, however, the 
complexity of the algorithm exponentially increases with respect to the 

25 maximum number of delayed symbols to be considered in the equalizer, 

resulting in very high computational complexity. For example, when BPSK 
is used as the modulation system, assuming that the storage capacity of one 



receiving channel is 10 symbols (the maximum number of delayed symbols to 
be considered in the equalizer being 10 symbols), the number of states 
becomes as large as 2'°=1024, exceeding practical computational complexity. 
The computational complexity could be reduced by using a linear equalizer 
for the iterative equalization. 

Fig. 2 illustrates a typical example of the linear equalizer. In the 
illustrated equalizer sampled received signals are applied to a 
cascade-connection of delay elements (of T-sec delay) 31 to obtain signals 
delayed by T-sec intervals in a sequential order, which are weighted by a 
variable weighting circuit 32 having set therein a coefficient by a coefficient 
adaptation circuit 34 and combined by an adder 33 to implement equalization. 
As the delay time x, the sampling period T or T/2 is usually chosen. 
Considering the equalizer as one filter, the linear equalizer achieves 
equalization by implementing a filter characteristic convenient for identifying 
the received signal. This filter basically implements an inverse characteristic 
of a channel distortion. 

When such a linear equalizer is used for the above-described iterative 
equalization, it is necessary to derive a soft decision value as the equalizer 
output. In the afore-mentioned literature by Daryl Reynolds and Xiaodong 
Wang, there is set forth a scheme in which, by deriving the soft decision value 
through the use of the linear equalizer, the iterative equalization can be 
achieved with less computational complexity than in the case of using the 
conventional algorithm such as MAP, SOVA or the like. In the literature a 
description is given of the application of the linear equalizer assuming the 
following discrete-time model as a channel model. A received signal sample 
r(n) is expressed as follows: 



r(n) = X t>(n - j)h(n; j) + cjv(n) , n=0, 1 , 
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(3) 



v(n) =[vo(n) vi(n) . . .VM-i(n)]'^ (4) 
h(n;j) =[ho(n;j) hj(n;j) . . .hM.i(n;j)]'' (5) 
In the above, M is the number of channel outputs and 

r(n)=[ro(n)ri(n). . .rM.i(n)]^ is a vector representing a received sample value in 
each channel at time n, where indicates a transposed matrix. For example, 
in the case of using a plurality of antennas, M corresponds to the number of 
antennas. B(n) represents a coded bit, and v(n) in Eq. (4) is a complex 
vector of an average 0, representing noise. In Eq. (5), h(n;j) is a vector 
representing the tap weighting coefficient of the channel, and J is the channel 
ISI length. Now, a received signal vector R(n), a noise vector V(n), a coded 
bit vector B(n) and a channel matrix B.Jn) are defined as follows: 
"r(n + J - 1)' 

R(n)= : (6) 

-MJxl 



V(n) = 



r(n) 
v(n + J - 1) 

v(n) 



(7) 



B(n)=[b(n+(J-l))b(n+(J.2)). . .b(n+l)b(n)b(n-l). . .b(n-(J-l))]T (8) 
H(n) = [h(n;0)-..h(n;J-l)] (9) 

'H(n) 0 0 

0 H(n) :* 

1 0 
^ - 0 H(n)J^^^,^_^^ 



Hn,(n) = 



h(n;0) h(n;l) 
0 h(n;0) h(n;l) 



h(n;J-l) 



h(n;J-l) 



0 



0 



0 



0 



(10) 



0 



0 



h(n;0) h(n;l) 
0 h(n;0) 



h(n;J-l) 



0 



0 



h(n;l) 



h(n;J-l) 



MJx(2J-l) 



The vectors R(n) and V(n) given by Eqs. (6) and (7) are column vectors 
of a matrix with MJ rows and one column, that is, MJ-dimensional column 
vectors. B(n) is a (2J-1 )-dimensional column vector, H(n) is a matrix with 
M rows and J columns, and UJn) is a matrix with MJ rows and (2J-1) 
columns. From the above, Eq. (3) becomes as follows: 

R(n) =Hnj(n)B(n) +aV(n) (1 1) 

Next, the likelihood of the code bit string is defined by the following 
equation using the a priori LLR value ?i2[b(k)] of the code bit string b(i). 
b'(k) =tanh[;t2[b(k)]/2]. n-(J-l) <k<n+(J-l) (12) 
Eq. (12) has the property that it rapidly approaches 1 or -1 upon 
deviation of a variable x of a hyperbolic function tanhx from 0, and the 
likelihood b'(k) represents an estimated value of b(k). Further, that 
estimated value vector of the code bit in the received signal vector R(n) which 
affects, as intersymbol interference, code bit b(n) at time n is given by the 
following equation. 

B'(n)=[b'(n+(J-l))b'(n+(J-2)). . .b'(n+l)Ob'(n-l). . .b'(n-(J-l))]^ (13) 
In Eq. (13), since B'(n) is an estimated value vector that is used to subtract 
from the received signal vector R(n) the intersymbol interference component 
for the code bit b(n) at time n, the element corresponding to b(n) of the vector 
B'(n) is set to zero. A difference vector Rc(n) by subtracting an intersymbol 
interference component replica from the received signal vector R(n) is 
defined by the following equation using the estimated intersymbol 
interference code vector B'(n). 
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Rc(n) = R(n)- Han)B'(n) = HUn)[B(n)-B'(n)]+av(n) (14) 
The difference vector Rc(n) is subjected to filtering expressed by the 
following equation through the use of a linear filter of tap coefficient G(n). 
Z(n) = G(n)«R,(n) (15) 
5 where " is a conjugate transposed matrix. The tap coefficient G(n) is so 
chosen as to minimize the mean square error between the code bit B(n) and 
the filter output Z(n) as expressed by the following equation. 
G(n)=arg min E{||b(n)-G(n)"Re(n)|p} 

=arg min G(n)"E{Rc(n)Rc(n)"}G(n)-2 G(n)"E{b(n)R,(n)} (16) 
10 where arg min E{} represents that one of all J-dimensional vectors G whose 
content of { } is minimum. 

Expansion of first and second terms in Eq. (16) gives 
E{Rc(n)R,(n)"}=H„(n)A(n)H^(n)"Wl (17) 
E{b(n)Rc(n)}=H^(n)ej (18) 
1 5 where I is a unit matrix and 
A(n)=Cov{B(n)-B'(n)} 

=diag[l-b'V(J-l))v..,l-b''(n+l),l,l-b''(n-lX 

...,l_b'V(J-l))] (19) 
In the above, Gov represents a covariance matrix and diag represents a 
20 diagonal matrix. Cj is a vector of a length 2J-1 , a J-th element being 1 and 
other elements being Os. By solving the normal equation with the gradient 
vector on the right-hand side of Eq. (16) to 0 and substituting Eqs. (17) and 
(18), the optimum tap coefficient G(n) is given as follows: 

G(n)=[H„,(n)A(n)Hn,(n)Wl]-'H„,(n)ej (20) 
25 Substitution of the thus derived tap coefficient G(n) into Eq. ( 1 5) gives 
Z(n)=e/H4n)«[H^(n)A(n)H„,(n)"+o'l]-' [R(n)-H^(n)B(n)] (2 1) 
In this case, since it is necessary that Hm(n) be a fiill column rank, the antenna 



-8- 



number M needs to be at least 2 or more. 

Assuming that the filter output can be approximated by a Gaussian 
distribution, Z(n) by Eq. (21) can be described as follows (see V. Poor and S. 
Verdu, "Probability of Error in MMSE Multiuser Detection," IEEE Trans. 
5 Information Theory., vol. IT-43 , No. 3 , pp. 85 8-87 1 , May 1 997). 
Z(n)=M,(n)b(n)+n(n) (22) 
where fi(n) is an equivalent amplitude of the output signal and r|(n) represents 
a Gaussian distribution of an average 0 and the variance v^(n). Therefore, 
^(n) and v(n) can be expressed as follows: 
10 Kn)=E{Z(n)b(n)} 

=[Hn,(n)«[H„,(n)A(n)H„,(n)"Wl]-'H^(n)]j,j (23) 
where indicates the element at the intersection of the J-th row and the J-th 
column of the matrix. 

=var {Z(n)}=n(n)-li'(n) (24) 
1 5 From the above, the extrinsic information derived by the linear equalizer can 
be deduced from the following equation. 

^p(Z(n)|b(n) = -l) l-n(n) ^^"^ 
To obtain the optimum tap coefficient G(n) with the above-described 
method, however, it is necessary that the inverse matrix calculation expressed 
20 by the following equation be conducted in Eq. (20) for each point in time; 
hence, much time is required for calculation. 

cD(n) = [H^(n)A(n)Hn,(n)" + ah]-' (26) 
With the conventional method for iterative equalization by the linear 
equalizer, the inverse matrix calculation needs to be conducted for each point 
25 in time to update the tap coeflficient~this gives rise to the problem of high 
computational complexity. 



SUMMARY OF THE INVENTION 

It is therefore an object of the present invention to provide an adaptive 
equalization method and apparatus that permit reduction of the computational 
complexity for deriving the tap coefficient of the linear equalizer for iterative 
equalization and hence allow shortening the time for computation. 

According to a first aspect of the present invention, a received signal 
sample value and a known signal are used to estimate the impulse response of 
the channel, and a received signal sample value sequence is used to adaptively 
calculate the tap coefficient of a linear filter. The tap coefficient is used to 
perform linear filtering of the received signal sample value sequence, and the 
result of linear filtering and the estimated value of the impulse response are 
used to derive a soft decision value. 

With this construction, in the case of using the linear equalizer for 
iterative equalization, the received signal sample value sequence is input 
thereto to adaptively derive the tap coefficient of the linear filter-this permits 
reduction of computational complexity as compared with the prior art method 
involving the inverse matrix calculation. 

According to a second aspect of the present invention, in second and 
subsequent processing by an adaptive iterative equalizer for the iteration of 
equalization and decoding, an estimated value vector of the code bit based on 
the likelihood of a soft decision value firom a decoder is subjected to linear 
filtering with an estimated impulse response value Hni(n) to generate a replica 
of the received signal. The replica is subtracted fi-om the received signal to 
create a difference signal Rc(n) having removed therefirom the intersymbol 
interference component, then the tap coefficient of the linear filter is 
calculated, for example, using the Matrix Inversion Lemma, fi-om the 
estimated impulse response value vector HlCh) approximated with the 
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intersymbol interference component regarded as zero, then the difference, 
then the difference signal Rc(n) is subjected to linear filtering using the tap 
coefficient, and the result of linear filtering and the estimated hnpulse 
response value vector are used to derive the soft decision value. 

Mth this construction, the tap coefficient of the linear filter in second 
and subsequent iterative equalization by the linear equalizer can be calculated, 
for example, using the Matrix Inversion Lemma, by which it is possible to 
achieve adaptive equalization with less computational complexity than in the 
case of the conventional that involves the inverse matrix calculation of Eq. 
(20). 

According to a third aspect of the present invention, a received signal 
sample value sequence and a knovm signal are used to estimate the impulse 
response of the channel, then a threshold value is set accordingly, and 
components of values smaller than the threshold value are removed from the 
estimated impulse response value vector, that is, the dimension of the 
estimated impulse response value vector is reduced. Further, components 
corresponding to those of the values smaller than the threshold value and 
hence removed from the estimated impulse response value vector are removed 
from the tap coefficient of a linear filter calculated from the received signal 
sample value sequence and the received signal vector to thereby reduce their 
dimension as well. The soft decision value is derived using the estimated 
impulse response value vector, the tap coefficient of the linear filter and the 
received signal vector of the thus reduced dimensions. 

The above scheme permits reduction of the computational complexity 
for deriving the soft decision value when the linear equalizer is used for 
iterative equalization. 

According to a fourth aspect of the present invention, the likelihood of 



the soft decision value from a decoder is calculated in second and subsequent 
processing by an adaptive equalizer for iteration of equalization and decoding, 
then the tap coefficient of the linear filter is calculated, for instance, using the 
Matrix Inversion Lemma from the estimated impulse response value vector 
reduced in dimension in the first round of iteration and stored in a memory, 
then the tap coefficient is used to perform linear filtering of a difference signal 
having subtracted therefrom an intersymbol interference component, and the 
result of linear filtering and the estimated impulse response value vector are 
used to derive the soft decision value. 

This scheme allows reduction of the computational complexity for 
deriving the soft decision value in the second and subsequent iterative 
equalization by the linear equalizer as well. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram of a transmitter-receiver of the iterative 
equalization scheme; 

Fig. 2 is a block diagram of a linear adaptive equalizer; 

Fig. 3 is a block diagram illustrating the fiinctional configuration for 
first equalization processing in a first embodiment of the adaptive equalizer 
according to the present invention; 

Fig. 4 is a block diagram of the fimctional configuration in second and 
subsequent rounds of equalization; 

Fig. 5 is a graph for explaining the effect of the present invention; 

Fig. 6 is a diagram depicting an example of a multipath channel 
environment; 

Fig, 7 is a block diagram illustrating the functional configuration for 
first equalization processing in a second embodiment of the adaptive equalizer 
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according to the present invention; 

Fig. 8 is a block diagram of the functional configuration in second and 
subsequent rounds of equalization; 

Fig. 9 is a graph showing an example of an impulse response 
5 contained in the received signal; 

Fig. 10 is a graph showing an example of selecting effective paths 
from the impulse response; 

Fig. 1 1 is a graph for explaining the effect of the present invention; 

and 

1 0 Fig. 1 2 is a graph for explaining a channel model used for computer 

simulations. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
EMBPDIMENT 1 

15 In the method set forth in the afore-mentioned literature that employs 

a linear equalizer for iterative equalization, the likelihood b'(k) of the code bit 
is calculated by Eq. (12) using the a priori LLR value A-2[b(k)] of the code bit 
string. In the first equalization processing, since the a priori LLR value 
^2[b(k)] in Eq. (12) is zero in the absence of the a priori information, the 
20 likelihood b'(k) is all zeros as shown below. 
b'(k) = [OO...O...f 
Therefore, the difference vector Rc(n) by Eq. (12) becomes as follows: 

Rc(n) =R(n) -Hn,(n)B'(n) =R(n) (27) 
and A(n) becomes a (2J-l)x(2J-l) whose diagonal components are all ones 
25 as shown below. 

A(n) = Cov{B(n)-B'(n)}=diag{l,l,...,l, ...,1,1} (28) 
Accordingly, the optimum tap coefficient G(n) of the linear filter by 
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Eq. (20) can be derived by the following calculation. 

G(n)= [H„,(n)H„,(nf +o'l]-^H„,(n)ej (29) 
From Eq. (29) the tap coefficient G(n) is dependent on the value HJn) 
corresponding to the impulse response of the channel. In the case v^^here 
temporal variations of the channel due to high transmission rate or the like, 
that is, where the duration of a received burst path is so short that the channel 
characteristic of the received signal does not change in that time interval, the 
tap coefficient G(n) is substantially constant. Accordingly, by precalculating 
the tap coefficient G(n) prior to the processing of the information sjrabol 
period and using the precalculated value for processing of each information 
symbol, the computational complexity can be reduced as compared with that 
in the conventional system that derives the tap coefficient G(n) for each 
information symbol. Moreover, since the above-mentioned temporal 
variations need not be taken into account as well in the calculation of the tap 
coefficient by Eq.(20), an adaptive algorithm can be used~this avoids the 
necessity for conducting the inverse matrix calculation to derive the tap 
coefficient G(n), permitting reduction of computational complexity. 

In view of the above, in the first embodiment of the present invention, 
the tap coefficient calculation is conducted by an adaptive algorithm to obtain 
G(n) in a known signal period or training signal period prior to the processing 
of the information symbol period, and the value thus obtamed is used to 
perform the processing in the information symbol period. This eliminates 
the needs for the inverse matrix calculation and the derivation of G(n) at each 
point in time of the information symbol, enabling the computational 
complexity to be made lower than in the conventional scheme. 

Figs. 3 and 4 depict examples of functional configurations for first 
processing and the second and subsequent processing for equalization by the 
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adaptive equalizer according to the first embodiment of the present invention. 
The illustrated adaptive equalizer corresponds to the equalizer 21 in Fig. 1. 
Fig. 3 shows the functional configuration for the first equalization processing, 
and Fig. 4 for the second and subsequent equalization processing. The 
broken line in Fig. 4 indicates the portion that is not used at the current point 
in time (in second and subsequent equalization processing) but used in the 
first equalization processing. 

The adaptive equalizer of the first embodiment comprises an impulse 
response estimating part 100, a linear filter part 1 10, a soft decision value 
calculating part 120, a switch 101, a likelihood calculating part 130, and a 
switch 160. The impulse response estimating part 100 is made up of a 
transversal filter 103, a subtracter 105, and an adaptive estimating part 106. 
The linear filter part 11 0 is composed of a transversal filter 11 1 , a tap 
coefficient calculating part 1 12, and a subtracter 113. 

A description will be given, with reference Fig. 3, of the first 
equalization processing in the iterative equalization. In this instance, the 
switch 101 is connected to a known signal, for example, training signal input 
terminal 102, and the switch 160 is coimected to an input terminal 104. In 
the impulse response estimating part 100, a training signal sequence b(n), 
provided as a known pattern fi-om the training signal input terminal 102, and a 
difference value (vector) Rc(n), obtained by subtracting an output replica 
Hm(n)b(n) of the transversal filter 103 by the subtracter 105 fi-om a received 
signal sample value sequence R(n) fi-om the input terminal 104 are used to 
estimate the impulse response Hn,(n) of the channel in a calculating part 106 A 
of the adaptive estimating part 1 06. This impulse response estimating 
calculation needs only to minimize the square of Rc(n) through the use of an 
adaptive algorithm utilized in an ordinary adaptive equalizer. The estimated 
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impulse response value Hm(n) having converged at the point of termination of 
the training signal sequence b(n) is stored in a storage part 106B. 

In the linear filter part 110, the received signal sample value sequence 
R(n) is input to the transversal filter 1 1 1 as a linear filter, by which it is 
subjected to linear filtering with the tap coefficient G(n) calculated in the tap 
coefficient calculating part 1 12, then the filtered output Z(n)=GH(n)"R(n) is 
subtracted by the subtractor 113 firom the training signal sequence b(n), and 
the difference is input to the tap coefficient calculating part 1 12 to update the 
tap coefficient G(n) of the linear filter. This updating calculation utilizes an 
adaptive algorithm. That is, in the first iteration of equalization, since 
R(n)=Rc(n) firom Eq. (27), the tap coefficient G(n) is so determined by the 
adaptive algorithm as to minimize ||b(n)-G(n)"R(n)|| firom Eq. 16. On this 
account, the tap coefficient G(n) can be calculated with less complexity than 
in the case of conducting the inverse matrix calculation by Eq. (20). The 
updating of the tap coefficient G(n) stops with the termination of the training 
signal sequence b(n), and the tap coefficient at this time is held unchanged. 

The received signal after termination of the training signal sequence 
b(n), lhat is, the information-symbol received signal sample sequence R(n), is 
subjected to linear filtering by the transversal filter 1 1 1 with the most recently 
updated tap coefficient G(n). The filtering output Z(n) and the ultimately 
estimated channel impulse response value Hm(n) are used to calculate, as the 
first equalization output, the soft decision value Xi[b(n)] in the soft decision 
value calculating part 120. In the first iteration of equalization, since the a 
priori information ^2[b(n)]=0, Ai=[b(n)]=Xi[b(n)]. The soft decision value 
A,i[b(n)] is fed to the de-interleaver 22 (see Fig. 1). 

Next, a description will be given first of the principle of equalization 
in the second and subsequent rounds of iterative equalization and then of 
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concrete equalization processing with reference to Fig. 4. In the second and 
subsequent rounds of equalization processing, the received signal R(n) used 
in the first iteration of processing is input again. On the other hand, when 
the likelihood b'(n) is calculated by Eq. (12) from the output of the decoder 23 
5 (Fig. 1) produced for the result of the first iteration equalization, to be more 
precise, from the soft decision value A,2[b(n)] of the code bit fed from the 
interleaver 25, the likelihood b'(n) is approximately +1 or -1 if Eb/No (Bit 
energy-noise ratio) is high. In other words, the estimated intersymbol 
interference code vector of Eq. (13) can be approximated as expressed by the 
10 following equation. 

B'(n)=[±l, ±1, ±1,0, ±1,..., ±1, ±lf (30) 
The use of thus approximated b'(n) permits approximation of Eq. (19) to the 
following equation. 

A(n)=Cov{B(n)-B'(n)}=diag{0,0,..., 0,1,0, ...,0,0} (31) 
1 5 Substituting Eq. (31) into Eq. (20), the tap coefficient G(n) of the linear filter 
can be obtained by conducting the following calculation. 

G(n)=[HL(n)HL«(n) WiY'H^ (32) 
where HL(n) is the vector of one column of M by J rows expressed by 

HL(n)=[h(n; J-l)...h(n; 0)]^=[ho(n; J-l)...hM-i(n; J-1) 
20 ho(n; J-2). . .hM-i(n; J-2). . .ho(n; 0). . .hM-i(n; 0)]"^. 

In this instance, the use of the Matrix Inversion Lemma in the 
calculation of the tap coefficient G(n) of the linear filter by Eq. (32) allows 
reduction of the computational complexity therefor. With such 
approximation of the estimated intersymbol interference code vector B(n) as 
25 expressed by Eq. (30), HL(n) becomes the vector of one column of M by J 
rows as mentioned above —this avoids the need for conducting the inverse 
matrix calculation, permitting substantial reduction of the computational 
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complexity as compared with that needed by the conventional scheme. Thus, 
according to the present invention, in the second and subsequent rounds of 
iterative equalization the computational complexity in the tap coefficient 
calculating part for deriving the tap coefficient G(n) of the linear filter can be 
reduced through utilization of the approximation by Eq. (3 1) and the Matrix 
Inversion Lemma. 

Fig. 4 illustrates in block form an example of the functional 
configuration in the second and subsequent rounds of equalization by the 
adaptive equalizer of the first embodiment. The change-over switch 101 is 
connected to the output of the likelihood calculating part 130, and the switch 
1 60 is connected to the output of the subtractor 105. The soft decision value 
A-2[b(n)] of each bit of the code sequence, which is the output from the 
decoder 23 (the output from the interleaver 25) in the previous round of 
iterative equation, is input to the likelihood calculating part 130 from a 
terminal 130, which conducts the calculation of Eq. (12), outputting the 
likelihood b'(n). The likelihood b'(n) is filter-processed by the transversal 
filter 103 using the estimated channel impulse response value Hm(n) stored in 
the storage part 106B, obtaining Hni(n)B'(n) in Eq. (14) as a replica of the 
intersymbol interference in the received signal R(n). The output replica of 
the transversal filter 103 is subtracted by the subfractor 105 from the received 
signal value sequence R(n) to derive the difference value Rc(n) of Eq. (14). 

The tap coefficient G(n) of the linear filter is calculated by Eq. (32) 
using, as the input therefor, only the estimated channel channel impulse 
response value HL(n) stored in the storage part 106B unlike in the scheme for 
deriving the tap coefficient of the linear filter in the first iteration of iteration. 

This calculation utilizes the afore-mentioned Matrix Inversion Lemma. 
That is, when this Lemma is expressed by A=B~'+CD"*C", where A and B 
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represent (M by M) positive definite matrixes, C an (M by N) matrix, D an (N 
by N) positive definite matrix, an inverse matrix of A is given by 

A"^=B-BC(D+C"BBC)'*C"B (33) 
Applying this Lemma to the calculation of the inverse matrix in Eq. (32), 
5 HL(n)HL"(n)+a2l=B-^+CD-'c", 

HL(n)HL"(n)=CD-'C", 

h(n) =C, 

10 HL"(n)=C" 

By using them to calculate Eq. (33), the inverse matrix calculation in Eq. (32) 
is solved. Incidentally, Eq. (33) also contains an inverse matrix calculation 
(D+C"bBC)~^ but this inverse matrix is also scalar, and hence it can similarly 
be calculated. 

15 The thus calculated tap coefficient value G(n) is used to perform 

linear filtering of the difference value Rc(n) by the transversal filter 111, and 
the filtering result Z(n) and the estimated channel impulse response value 
HL(n) stored in the storage part 106B are used to calculate Eq. (25) in the soft 
decision value calculating part 120 to derive the soft decision value Ai[b(n)]. 
20 The equivalent amplitude \i(n) in Eq. (25) is given by the following equation 
by substituting Eq. (31) for A(n) in Fig. (23). 

Kn) = [HL(n)«[HL(n)HL(n)« +<^lfHi^(n)] 
=HL(n)"G(n) (34) 
Incidentally, in the first iteration of equalization the tap coefficient calculating 
25 part 112 of the linear filter part 110 may also use the Matrix Inversion 
Lemma to calculate the tap coefficient G(n). 

As described above, this embodiment permits reduction of the 
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computational complexity for deriving the tap coefficient of the linear 
equalizer that is used for iterative equalization. Fig. 5 shows an error rate 
characteristic obtained by computer simulations, for example, when an RLS 
(Recursive Least Squares) algorithm of excellent convergence characteristic is 
5 used as an adaptive algorithm for deriving the tap coefficient of a linear 
equalizer used for iterative equalization. 

In Fig. 5, the ordinate represents BER (Bit Error Rate) and the 
abscissa Eb/No (bit energy-noise ratio). The modulation system used was 
BPSK (two-phase shift), the information transmission rate was 12 Mbps, and 

10 the maximum Doppler frequency indicating the fading variation rate was 

1000 Hz. The frame configuration used is a unique word of a known pattern 
for first 128 symbols, then the channel impulse response Hm(n) was estimated, 
and 128 symbols were added as information symbols to the end of the training 
symbols. The RLS algorithm was used for the estimation of the impulse 

1 5 response. The channel was set as a five-wave Rayleigh fading channel of 
the same level, and the receiving side was to perform two-branch diversity 
and iterative equalization by a five-tap linear filter. Further, SOVA 
(Soft-Output Viterbi Algorithm) was used for the decoding in the iterative 
equalization, and the number of iterations was two. 

20 In Fig. 5, the curves AO and Al show tap coefficients of the linear 

filter obtained by conducting the inverse matrix calculation at respective point 
in time in the cases where the number of iterations were one and two, 
respectively. The curves BO and B 1 show the tap coefficients G(n) derived 
by the scheme according to the first embodiment of the present invention, the 

25 curve BO indicating the tap coefficient values by using the adaptive algorithm 
and the curve Bl the tap coefficient values by using the Matrix Inversion 
Lemma when the number of iterations was two. In the case of the present 
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invention, degradation of the bit error ratio is only on the order of 0.1 dB as 
compared with the ratio in the case of conducting the inverse matrix 
calculation at each point, it can be regarded as about the same as in the latter 
case. The present invention does not need to conduct the inverse matrix 
5 calculation expressed by Eq. (26), and hence it permits substantial reduction 
of the computational complexity for deriving the tap coefficient of the linear 
equalizer that is used for iterative equalization. Incidentally, it was 
confirmed that the same goes for the cases where the known pattern is 
128-symbol and the number of information symbols is 256 and where the 

3 10 known pattern is 1 28-symbol and the number of information symbols is 5 1 2. 

I This means that when the duration of a burst path (the frame length) is short 

P about the same as these information symbol durations, impulse response 

Q 

; variations of the channel can be ignored, and the bit error rate is 

i approximately equal to that in the case of conducting the inverse matrix 

J 15 calculation at each symbol time, and with less computational besides. 

J Letting M and J represent the number of antennas (channels) and the length of 

^ intersymbol interference (the number of taps of the linear filter 111), 

respectively, the computational complexity is around 0((MJ)^) in the case of 

conducting the inverse matrix calculation at each symbol time, but in the 
20 above-described embodiment of the present invention the computational 

complexity is approximately 0((MJ)^), which is substantially lower than in the 

above case. 

While the examples of Figs. 3 and 4 are shown to include two 
transversal filters 103 and 111 with a view to facilitating a better 
25 understanding of the function configuration, the number of filters need not 
always be two but may also be one. Moreover, the tap coefficient 
calculating part 112 and the adaptive estimation part 106 are shown as 
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different blocks, but they can be combined into one unit. 

As described previously, the one mode of the first embodiment utilizes 

the fact that the impulse response of the channel does not so vary over a short 

range of time, and hence it is characterized in that the training signal at the 
5 beginning of the received burst path (frame) or the known signal is used to 

obtain the estimated impulse response value Hni(n), which is used to the 

subsequent processing. 

Further, according to the other mode of Hiis embodiment, in the 

second iterations of iteration, only a J-th element of a diagonal matrix is set to 
□ 10 1 which is a covariance matrix of the difference between the code bit vector 
I B(n) and the estimated intersjnnbol interference code vector B'(n) as 

^ expressed by Eq. (3 1), and the other diagonal elements are approximated to 

'1 zeros, by which the calculation for the tap coefficient G(n) becomes a vector 

^ calculation as indicated by Eq. (32)~this avoids the necessity of conducting 

^ 15 the inverse matrix calculation and hence allows reduction of computational 
y complexity. In the first processing the known signal (training signal) need 

" not always be used, but instead the output from a discriminating part 150, 

which discriminates the output from the transversal filter 103 between two 

value in binary, may be used in place of the signal from the output of the 
20 switch 101 as indicated by the broken lines in Fig. 3. The estimated impulse 

response value Hin(n) ultimately obtained by the first processing is used in the 

second and subsequent processing. 

EMBODIMENT 2 

25 In the first embodiment described above with reference to Figs. 3 and 

4, all those of path components so-called multipath received signal R(n) 
arriving at the input terminal 104 via different paths which are withm the 
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range of the maximum number J of delayed symbols to be considered are 
equalized as will be seen from Eqs. (5), (6), (8) to (11), (15) and (25). 
According to the channel environment of the channel, however, path 
components of relatively large received power may sometimes arrive apart in 
time from other path components as indicated by multipath components dl, 
d2, . . . in Fig. 6. In such a channel environment, too, the first embodiment 
performs equalization of all the received path components concerned. 

That is, in the first embodiment, path components within the range of 
the maximum number J of delayed symbols, even if very small, in extreme 
cases, zero in received power, are also subjected to equalization, and 
consequently, much time is needed for calculation. As a solution to this 
problem, this embodiment is adapted to omit particular ones of the 
dimensions corresponding to multipath components of negligibly small 
received power to reduce the sizes of vectors and matrixes to be handled, 
thereby reducing the computational complexity and hence shortening the 
calculating time. 

Calculating the likelihood b'(n) of the soft decision value A,2[b(n)] for 
the code bit string b(n) that is the decoder output for the result Ai[b(n)] of the 
first iteration of equalization, by Eq. (12) as in the second iterations of 
equalization in the first embodiment, the calculated value is virtually +1 or -1. 
Through utilization of this fact for the approximation of the estimated 
intersymbol interference code vector B'(n) as expressed by Eq. (30), the tap 
coefficient G(n) of the linear filter 111 can be calculated using the Matrix 
Inversion Lemma as follows: 

G(n)= [HL(n)HL(n)" + a'l]-'HL(n) 

= HL(n)/(a^ + HL(n)"HL(n)) (35) 
where HL(n) is a column vector of a M by J length as shovra 
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below. 

HlCh) = [ho(n; J-l)...hM-i(n; J-1) ho(n; J-2)...hM-i(n; J-2) ... 
...ho(n;0)...hM.i(n;0)]^ (36) 
By the approximation of A(n) in Eq. (23) as indicated by Eq. (32), the 
5 equivalent ampUtude [i{n) of the output signal is also calculated as expressed 
by (34). 

On the other hand, Z(n) is given by 

Z(n)=G(n)"Rc(n) (37) 
as is the case with Eq. (15). By using in Eq. (25) the linear filter output Z(n) 
10 obtained by Eq. (37), the extrinsic information value Ai[b(n)] derived by the 
linear equalizer can be used as the results of second and subsequent 
equalizations. 

Now, assume that path components j=l to j=J-2 within the range of the 
maximum number J of delayed symbols to be considered in the equalizer have 
1 5 not arrived but that path components j=0 and j=J- 1 have arrived. If the 
impulse response of the channel can accurately be estimated, other path 
components than effective ones are such that h(n;l),h(n;2),...,h(n;J-2)=0. 
Accordingly, it follows from Eq. (35) that 

HL(n)=[ho(n;J-l)...hM-i(n;J-l) 0...0 0...0 ho(n;0) 

20 ...hM-i(n;0)f (38) 

Since the tap coefiScient G(n) of the linear filter can be calculated 
from HL(n) as shown by Eq. (35), it follows that 

G(n)= [go(n;J-l). . .gM-i(n;J-l) 0. . .0 go(n;0). . .gM-i(n;0)]^ (39) 

The derivation of the soft decision value Ai[b(n)] by the linear 
25 equalizer by Eq. (25) requires the calculation of the equivalent amplitude |j,(n) 
of the output signal by Eq. (36) and the linear filter output Z(n) by Eq. (37). 
The equivalent amplitude |i(n) and the linear filter output Z(n) are calculated 
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as shown below. 

Kn) = HL(n)"G(n) 

= ho(n; J-l)xgo(n;J-l)+...+ hM-i(n; J-l)xgM-i(n;J-l)+...+OxO+... 
...+0x0+...+ ho(n;0)xgo(n;0)+...+ hM-i(n;0)xgM.i(n;0) 

(40) 

Here, in the case where zero-elements of HL(n)" and G(n) are omitted, 
that is, where the dimensions of HL(n)" and G(n) are lowered as shown below, 
too, the value of fj,(n) naturally remains unchanged. 

HL'(nf = [ho(n;J-l). . .hM-i(n;J-l) ho(n;0). . .hM-,(n;0)] (41) 
G'(n) = [go(n;J-l)...gM-i(n;J-l) go(n;0)...gM-i(n;0)f (42) 
^(n) = HL'(n)"G'(n) 

= ho(n;J-l)xgo(n;J-l)+. . .+ hM-i(n;J-l)xgM.,(n;J-l) 

+ho(n;0)xgo(n;0)+. . .+hM-i(n;0)xgM.i(n;0) (43) 
And, Z(n) becomes as follows: 
Z(n) = G(n)"Rc(n) 

- go(n; J-l)xro(n+J-l)+...+ gM-i(n; J-l)xrM-i(n+J-l)+... 

+Oxro(n+J-2)+. . .+0xrM-i(n+J-2)+ +Oxro(n+l)+. . .+0x 

rM-i(n+l)+go(n;0)xro(n)+. . .+gM-i(n;0)xrM-i(n) (44) 
In Eq. (44), since G(n)" has 0-elements, the values of Rc(n) 
corresponding to the 0-elements are unnecessary in the derivation of Z(n) as 
shown below. 

Rc'(n) = [ro(n+J-l). . .rM-i(n+J-l)ro(n). . .rM-i(n)] (45) 
Z'(n) = G'(n)»Rc'(n) 

= go(n;J-l)xro(n+J-l)+...+ gM-i(n; J-l)xrM-i(n+J-l)+... 
. . .+ go(n;0)xro(n)+. . .+ gM.i(n;0)xrM.i(n) (46) 
Since the dimension of Rc(n) can be lowered as described above, the 
dimension of the received signal vector R(n) can also be lowered by Eq. (14). 



-25- 

Further, since the dimension of the intersymbol interference component 
replica Hn,(n)B'(n) can also be reduced, the dimension of the channel matrix 
Hin(n) can be reduced as well. 

From the above, it can be seen that calculations for the path 
5 components having arrived can be dispensed with in the calculation for 
deriving the extrinsic information by the linear equalizer. That is, in the 
estimated impulse response value of the channel, the tap coefficient G(n) of 
the linear filter, the received signal vector R(n) and Rc(n) obtained by 
subtracting the intersymbol interference component from the received signal 
10 vector, those elements of vectors and matrixes corresponding to path 

components having not arrived can be removed. Accordingly, it is possible 
to decrease the number of dimensions of vector components and the number 
of elements of matrixes necessary for the calculation involved in the 
derivation of the soft decision value of the linear equalizer. 
1 5 However, the estimated impulse response value does not become zero 

due to an error in the estimation of the impulse response of the channel. 
Hence, a certain threshold value hth is set for comparison and discrimination 
of the estimated impulse response value as described below. For example, if 
h(n;l),h(n;2),...h(n;J-2)<hth, then it is regarded as 
20 h(n; 1 ),h(n;2), . . . h(n; J-2)=hth. Letting L represent the number of effective 
received path components larger than the threshold value hth, the dimension of 
the estimated impulse response vector HL(n), which indicates all estimated 
impulse response values within the range of the maximum number J of 
delayed symbols, can be lowered from an M* J-th dimensional vector to an 
25 M*L-th dimensional vector H'lCu) as expressed below. 

HL(n) =[ho(n; J-l)...hM-i(n; J-1) 0 ... 0 ho(n;0)...hM-i(n;0)f (47) 
HL'(nHho(n; J-1). . .hM-i(n; J-l)ho(n;0). . .hM-i(n;0)]'^ (48) 
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Further, the number of row elements of the channel matrix Hni(n) can 
be lowered from 2J-1 to 2L-1. 



Hm(n) = 



h(n;0) h(n;l) 



h(n;J-l) 



h(n;0) h(n;l) 



h(n;0) 
0 



h(n;l) 
h(n;0) 



0 

h(n;J-l) 0 

h(n;J-l) 



h(n;l) 



0 



0 

h(n;J-l) 



MJx(2J-l) 



(49) 



HmXn) 
If 



_rh(n;0) 0 ••• h(n;J-l) 0 ••• 0 ~| 
~\_ 0 — 0 h(n;0) 0 ••• h(n;J-l)J 



1mLx(2J-1) 



(50) 



h(n; l),h(n;2),. . .h(n;J-2)<hth, 



then 



go(n; 1). . .gM-i(n; l),go(n;2). . .gM-i(n;2),. . .,go(n;J-2). . .gM-i(n;J-2)= 0 
is set; the tap coefficient vector G(n) of the linear filter in the second 
iterations of iteration is given by the following equation. 

G'(n)= [go(n;J-l)...gM.i(n;J-l) go(n;0)...gM-i(n;0)]^ (51) 
This is apparent from Eq. (35) in the calculation of HL'(n)"HL'(n). Moreover, 
the tap coefficient vector G(n) can also be lowered from the M* J-th 
dimensional vector to the M*L-th dimensional vector G'(n). 

Similarly, the received signal vector R(n) and the difference vector 
Rc(n) can also be lowered from the M* J-th to the M*L-th dimension. 

The equivalent amplitude |i(n) of the output signal can be calculated 
by Eq. (43) from HlXh) and G'(n), and the linear filter output Z'(n) can be 
calculated by Eq. (46) from G'(n)*^ and Rc'(n). Accordingly, they can each 
be derived by calculations of vectors of low dimensionality. 

Furthermore, it is possible to lower the dimension of the number of 
column elements at positions of path components regarded as zeros in the 
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channel matrix Hm'(n) of Eq. (50) generated from the estimated impulse 
response value as expressed by the following equation, and the dimension of 
the channel matrix for deriving the soft decision value by Eq. (25) can also be 
reduced. In an example in which only paths of j=0 and j=J are set as 
5 effective paths, the number of column elements can be decreased from (2J-1) 
to (2L-1). 

H^.(,) rh(n;0) h(n;J-l) 0 1 
™^ ^ L 0 h(n,0) h(n;J-l)Jj^,(2L-l) 

Figs. 7 and 8 illustrate in block form examples of functional 
configurations in the first iteration and second and subsequent rounds of 

10 equalization. In this embodiment, an identifying part 140 and a switch 170 
are additionally provided in the functional configuration of the first 
embodiment shown in Figs. 3 and 4, and the tap coefficient calculating part 
1 12 is composed of a calculating part 1 12 A, a storage part 11 2B and an 
identifying part 1 12C. The adaptive estimation part 106 is also added with 

15 an identifying part 106C. 

In the first place, the switch 101 is connected to the terminal 102, and 
the switches 160 and 170 are both cormected to the terminal 104. In the first 
iteration of equalization the processes for obtaining the estimated impulse 
response value Hni(n) and the tap coefficient G(n) are the same as in the first 

20 embodiment. Further, in Fig. 7, effective paths are identified by the 

identifying part 106C based on the estimated impulse response value Hn,(n) 
having converged at the end of the training signal. Assume, for example, 
that estimated impulse response values are obtained at respective symbol 
delay times within the range of the maximum number J of delayed symbols as 

25 shown in Fig. 9. The effective path identifying method is, for instance, to 
select paths above a certain threshold value hth as depicted in Fig. 10. The 
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threshold value hth is for instance, the average power of all received paths, and 
received paths above the average power are selected. The effective paths 
thus identified are each represented, for example, by the number of delayed 
symbols and stored in the storage part 1 06B together with the corresponding 
5 tap coefiRcient value h. This data may be stored as data of the vector Hl'(N) 
given by Eq. (48) or as data of the matrix Hm'(n) given by Eq. (48); and it can 
easily be obtained from the one to the other and vice versa. Further, only 
those of the tap coefficient vectors G(n) of the linear filter 111 which 
correspond to the effective paths identified by the identifying part 106C are 
u 10 each stored as the tap coefficient vector G'(n) in the storage part 1 12B. 

1=^ Moreover, the switch 160 is switched to the identifying part 140, in 

"■'4 

in which only an information symbol received signal R'(n) corresponding to the 

M' effective path in the signals received after the training signal, that is, in the 

I* information symbol received signals R(n), is selected, and the thus selected 

I i 15 signal is subjected to linear filtering by the transversal filter 111 using its tap 

iy 

13 coefficient G'(n) stored in the storage part 1 12B most recently. The linear 

i2 

filtering output Z'(n) and the most recently estimated impulse response value 
H'L(n) are used to calculate the soft decision value Ai[b(n)] in the soft 
decision value calculating part 120. 

20 In the second and subsequent rounds of equalization, the switch 101 is 

connected to the likelihood calculating part 130, the switch 160 to the 
identifying part 140 and the switch 170 to the output of the subtractor 105 as 
depicted in Fig. 8. The likelihood b'(n) representing the estimated 
intersymbol interference code, provided fi-om the likelihood calculating part 

25 130, is subjected to filtering by the transversal filter 103 using the estimated 
impulse response value HL(n) stored in the storage part 106B, and the output 
replica HL(n)b'(n) is subtracted by the subtractor 105 fi-om the received signal 
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sample value sequence R'(n) to obtain the difference value Rc'(n)- For only 
the path identified in the identifying part 106C of the impulse response 
estimating part 100, the difference value Rc'(n) is identified by the identifying 
part 140, and the identification result is input to the linear filter 111. The tap 
5 coefficient G'(n) of the linear filter 1 1 1 is derived from only the estimated 
impulse response value Hl'Ch) stored in the storage part 106B, by Eq. (35) 
through utilization of the Matrix Inversion Lemma, and the tap coefficient 
G'(n) is stored in the storage part 1 12B. The thus stored estimated impulse 
response value HL'(n), the tap coefficient G'(n) of the linear filter 111 and the 

10 identified difference value Rc'(n) are used to calculate the soft decision value 
Ai[b(n)] in the soft decision value calculating part 120. 

As described above, the second embodiment permits reduction of the 
computational complexity for deriving the soft decision value of the linear 
equalizer. Fig. 1 1 shows the results of computer simulations of the bit error 

15 rate characteristic in the case of using the RLS (Recursive Least Squares) 
algorithm of excellent convergence characteristic for deriving the tap 
coefficient of the linear equalizer in the first iteration of iterative equalization 
and the Matrix Inversion Lemma in the second and subsequent iterations of 
iterative equalization. 

20 The modulation system used is BPSK; the information transmission 

rate is 10 Mbps; the maximum Doppler frequency indicating the fading 
variation is 1000 Hz; and the frame configuration is one that a training signal 
of a known pattern of 63 symbols is added as an information signal to 256 
symbols. A channel, in which only path components j=0 and j=4 as 

25 indicated by the solid lines in Fig. 12 had arrived, was assumed as a channel 
model. The receiving side performs two-branch diversity and uses SOVA 
(Soft-Output Viterbi Algorithm) for decoding. On the assumption that the 
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channel was ideally estimated, the number of iterations was set to three. 

In Fig. 11, the curves Al, A2 and A3 show the results obtained in the 
cases where: equalization was performed once and iterated twice and three 
times, respectively; the tap number of the linear equalizer was considered for 
5 0 to 4T delayed paths; and the entire information of five taps was used as the 
tap coefficient of the linear equalizer. The curves B 1 , B2 and B3 show the 
results obtained in the cases where: equalization was similarly performed 
once and iterated twice and three times, respectively; and two-tap (only 0 and 
4T delayed paths) information was used as the tap number of tiie linear 
10 equalizer. 

As is evident from Fig. 1 1 , a comparison of the results obtained with 
the use of five-tap and two-tap information as the tap number of the linear 

equalizer indicates that two and three iterations of equalization provide about 
the same error rate characteristic with low computational complexity. 

15 

EFFECT OF THE INVENTION 

As described above, the first embodiment of the present invention 
enables the linear equalizer for iterative equalization enables to reduce the 
computational complexity for calculating the tap coefficient of the linear filter 

20 and hence shorten the time for the calculation. 

According to the second embodiment of the invention, in a ch^mel 
environment in which path components of relatively large power have arrived 
apart in time from other path components, it is possible to adaptively reduce 
the computational complexity for calculating the extrinsic information by the 

25 linear equalizer that performs iterative equalization. 



